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Abstract 

We explain how to exploit systematically the structure of nilpotent orbits to ob- 
tain a solvable system of equations describing extremal solutions of (super-)gravity 
theories, i.e. systems that can be solved in a linear way. We present the procedure 
in the case of the STU model, where we show that all extremal solutions with a flat 
three-dimensional base are fully described with the help of three different nilpotent 
orbits: the BPS, the almost-BPS and the composite non-BPS. The latter describes 
a new class of solutions for which the orientation of half of the constituent branes 
have been inverted with respect to the BPS one, such that all the centres are in- 
trinsically non-BPS, and interact with each others. We finally recover explicitly 
the ensemble of the almost-BPS solutions in our formalism and present an explicit 
two-centre solution of the new class. 



1 Introduction 



One of the great success of string theory has been to be able to provide a statistical inter- 
pretation to the Bekenstein-Hawking entropy of BPS black holes through the counting 
of D-branes in the weakly coupled regime [1]. The validity of the computation is ensured 
by supersymmetry, but it has nevertheless been proposed that this property could gen- 
eralise to non-BPS extremal black holes [2, 3, 4]. The classification of supersymmetric 
composite black hole solutions has permitted to understand the mismatch between the 
enumeration of spherically symmetric BPS black holes and the counting of BPS states 
within weakly coupled string theory [5]. Understanding the space of states associated to 
extremal black holes therefore clearly requires to have a global understanding of compos- 
ite extremal black hole solutions. In the recent years, a lot of techniques developed in the 
context of supersymmetric solutions have been adapted to the case of non-BPS extremal 
ones (see [6]-[25] for part of the literature). Indeed, it has been understood that a lot 
of features are in fact intrinsically more related to extremality than supersymmetry. In 
particular, the underlying system of equations is then solvable, which means that it can 
be solved in a linear fashion [12]. Such a linear structure is the key point for solving 
explicitly the equations, and this remark permitted to construct a lot of new non-BPS 
solutions [10, 11, 12]. However, it has up to now remain unclear where such a solvable 
system was coming from, and thus how to generalise the approach to find other solvable 
systems. The aim of the present paper is to address this issue. 

Since more than twenty years now, it is well-known that the stationary solutions of 
supergravity theories coupled to abelian vector fields and scalar fields parametrizing a 
symmetric space G 4 /K 4 are described by a non-linear sigma model coupled to Euclidean 
three-dimensional gravity, which scalar field V parametrize a symmetric space G/K* 
[26]. In the cases of interest (such has Kaluza-Klein supersymmetric theories without 
hyper-multiplets), G is a simple group and K* is a non-compact real form of its max- 
imal compact subgroup. The scalar momentum P is defined as the component of the 
Maurer-Cartan form V~ l dV in the coset component p = g t* of the Lie algebra. Solu- 
tions describing spherically symmetric black holes are then determined by the associated 
Noether charge in the Lie algebra g of G, and can therefore be classified in terms of G- 
orbits [26, 27]. In the extremal limit, the Noether charge is nilpotent and the spherically 
symmetric extremal black hole solutions are classified in terms of a sub-class of nilpotent 
orbits of G in g [28, 29]. It is shown that P then lies in the same nilpotent orbit, and 
more precisely in its intersection with the coset component p. It has been exhibited that 
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the first order system of differential equations defining the composite BPS black hole 
solutions can itself be derived from the property that the scalar momentum P lies in the 
nilpotent orbit associated to spherically symmetric BPS black holes [30]. We show in this 
paper that this property extends to all extremal composite black hole solutions with a 
flat three-dimensional base metric, 1 without assuming any spherical symmetry, nor even 
axisymmetry. In other words, they are solutions of a solvable system of differential equa- 
tions, which is defined by the nilpotent orbit in which P is constrained to lie in. We will 
be able to obtain all the possible solvable systems of a given theory from the study of its 
nilpotent orbits, and also to explore in a systematic way the structure of the solutions 
space. In the present paper, for the sake of clarity, we only study the case of the STU 
model, i.e. H — 2 supergravity in four- dimensions coupled to three vector multiplets 
with a cubic prepotential. The ideas explained here are nevertheless more general, and 
we intend to present the more general case of M = 8 supergravity in an upcoming paper. 
The generalisation to any Af = 2 supergravity coupled to vector multiplets associated 
to a very special Kahler geometry is straightforward. For the STU model, we show that 
all solutions can be described with only three different nilpotent orbits, leading to three 
inequivalent solvable system of differential equations. The first one is the Denef system 
of equations for BPS solutions [5, 31] and the second is the almost-BPS system [9]. We 
obtain the latter from the nilpotent orbit approach, and show that it allows us to recover 
all the previously known solutions of that class. The third one is new, and we call it the 
composite non-BPS system. In terms of type IIA supergravity, one can understand the 
almost-BPS system as a system of floating D-branes for which the constituent branes of 
one type have an inverted orientation, e.g. D6-D4-D2-D0. Similarly, the composite non- 
BPS system is associated to a system of D-branes for which the constituent branes of two 
types have their orientation inverted, e.g. D6-D4-D2-D0. We present here the compos- 
ite non-BPS system of equations, but, again for the sake of simplicity, only present one 
two-centres solution of this class. This example will be enough to exhibit key properties 
of the new class, in particular the fact that the centres are interacting, and to show that 
more general solutions exist with an arbitrary number of centres. 

An important point in understanding extremal solutions is the following. A generic 
single-centre solution of both the almost-BPS and the composite non-BPS system will be 
singular, and thus not physical. In order to be regular, the momentum P must fall down 
at any black hole horizon, either in the BPS or the so-called physical non-BPS orbit, 
which are associated to extremal single-centre solutions. However, regularity allows P to 

1 Excluding the extremal Kerr solutions and composite generalisations thereof. 
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lie in more general higher rank nilpotent orbits on a generic point of the base space. This 
is schematically depicted in Figure 1. It is crucial to understand that this is precisely the 
possibility for P to lie in higher rank orbits, where a single-centre solution would not be 
regular, that permits to obtain interacting non-BPS multi-centres solutions. Almost-BPS 
solutions can describe multi-centre interacting solutions [10, 11], but no interaction can 
take place between two non-BPS centres. It is of importance that the composite non-BPS 
multi-centre black holes that we exhibit in the present paper have interactions between 
genuinely non-BPS centres. This is an important step toward the understanding of the 
structure of the moduli space of non-BPS multi-centre solutions. 

Note also that the existence of such solutions is not in contradiction with the previ- 
ously obtained conclusion that there were no composite non-BPS solutions associated to 
the first order system of differential equations describing non-BPS single-centre solutions 
[19]. Indeed we find that the most general regular solution of the 'physical' non-BPS 
system (associated to the 'physical' non-BPS nilpotent orbit) are single-centre solutions 
with a possible bounded angular momentum, whereas the composite non-BPS solutions 
only exist in a more general system, the composite non-BPS system. One main differ- 
ence with the BPS solutions is, whereas for the latter the ADM mass and the flow of 
the scalar fields in the asymptotic region of a supersymmetric space-time are entirely 
determined by the asymptotic central charges Z(q,p), Zi(q,p), they also depend on the 
specific structure of the interior space-time for a composite non-BPS solution. This can 
be interpreted within the attractor mechanism by a lifting of the flat directions normally 
associated to single-centre non-BPS black holes. In the presence of interactions, the flow 
of the scalar fields in the asymptotic region is indeed governed by the 'auxiliary field' de- 
pendent 'fake superpotential' described in [21], such that instead of being determined by 
extremizing the 'fake superpotential', the 'auxiliary fields' associated to the flat directions 
are determined in function of the specific structure of the interior space-time. 

The paper is organised as follows. In section 2 we present, from a purely group the- 
oretical point of view, how nilpotent orbits generically define solvable systems, and then 
focus on the case of the STU model where we obtain three solvable systems that encode 
all extremal solutions of the STU model with a flat base. We restrict the arguments 
based on group theory to this section, and the reader not interested in the mathematical 
details can consistently skip this part, its results being recalled when needed in the rest 
of the paper. Section 3 is devoted to set up the theory and the conventions, in particular 
to relate the three-dimensional and four-dimensional quantities. In this section we also 
provide the expression of the ADM mass as a function of the asymptotic central charges 
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Physical non-BPS 




Figure 1: This figure depicts a simplified version of the Hasse diagram of the nilpotent orbits 
of SO (4, 4), where the orbits have been identified when E 8 rg\ conjugate in e§(8)- It describes 
the topology of the space of nilpotent elements in so (4, 4), such that two nilpotent orbits are 
connected by a line if the left orbit lies in the topological closure of the right one, see [32] 
for details. The lower rank orbits, on the left, corresponds to 1/2, 1/4 and 1/8-BPS solutions 
of M = 8 supergravity with a null entropy. Then, two orbits have the same rank, they are 
the 1/8-BPS and the physical non-BPS orbit associated to the regular extremal single-centre 
solutions. On the right side of the regularity line (dotted red line), the rank of the orbits are too 
high to describe regular single-centre solutions, but still allows for multi-centre configurations. 
In this paper, we study in details the three blue shaded orbits of this diagram. All the extremal 
solutions can be obtained from them as subcases. 

and phase parameters determined by the specific configuration of the interior solution. In 
section 4, we choose a duality frame to show how, from our point of view, we reobtain all 
known almost-BPS solutions, while section 5 is devoted to the new composite non-BPS 
system. In this two cases, we also show that the Ehlers rotation, which is part of the 
three-dimensional duality group, but not of the four-dimensional one, does not give any 
new solutions. Finally, we construct and analyse a complete two-centre solution of the 
composite non-BPS class in section 6. Some details of the conventions are relegated to 
an appendix. 



4 



2 Solvable subalgebras and equations of motion 
2.1 General solvable subalgebra 

One interesting feature of the almost-BPS system of differential equations [9, 10], is that, 
although not free, it admits a graded structure that allows to obtain each function as the 
solution of a free Laplace equation with a source term determined as a non-linear function 
of the lower grade functions. This renders the system exactly solvable, i.e. it amounts 
to recursively solve linear Laplace equations with known non-linear source terms. We 
explain in this section how this behaviour is in fact generic for extremal systems. In 
terms of the non-linear sigma model over G/K* coupled to Euclidean gravity in three 
dimensions, this graded structure originates naturally from the graded structure of a 
solvable subalgebra. 

Consider a solvable subalgebra n C Q, which satisfies by definition that 

3n6N, ad n ra n={0}, (2.1) 

where 

ad n n=[n,n] (2.2) 
as a set, and the power n defines the number of commutators, n admits a grading 

n 

n = J>«, (2.3) 
P =i 

such that 

n (p) ad t f- 1 n\ad,fn . (2.4) 
We will consider that this grading is consistent with the involution defining the subalgebra 

r 

g = r©p (2.5) 

associated to K*, such that each component decomposes accordingly into 

n (p) { « pM . (2.6) 

Now, consider the Ansatz for the scalar field 

V = exp(-L) (2.7) 
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defined such that L is a function of the three-dimensional base M 3 valued in n fl p 



L : M 3 -> n n p . (2.8) 

L is in n because we assume it to be nilpotent, and it can be chosen in p by fixing the 
coset representative V in G/K* to be odd with respect to the involution defining K* ? 
One straightforwardly obtains that the components of the Maurer-Cartan form 

V~ 1 dV = P + B, Pep, Bet (2.9) 



arc 



P = -V ; ; 1 r . ad T 2k dL , B = -Y , , 1 x, ad f 2fc+1 dL , (2.10) 

k>0 K ' k>0 K ' 

and lie respectively in 

n n 

P ennp = J^p (p) , Bennf = ^f w . (2.11) 

p=i p=i 

Note that both sums are finite since ad L n dL = 0. By property of the solvable algebra 

Tr P^Py = . (2.12) 

Equations of motion. The equations of motion of (super-)gravity theories coupled to 
abelian vector fields in four dimensions, reduce for stationary solutions to the following 
equations on the three-dimensional Riemannian space M 3 

TtP^P v = R^ , (2.13) 
d*P+[B,*P] = . (2.14) 

So (2.12) translates into = 0, and the base three-dimensional metric is flat 7^^ = 5^. 
Therefore the base manifold is simply R 3 , although it will be more precise to define it as 
the punctured R 3 

M 3 = R 3 \{x. 4 |Li} , (2.15) 



2 Note that this choice of parametrization of the coset space G/K* is not globally defined, because 
G/K* is not topologically trivial as opposed to G/K. Nevertheless it is well defined on a dense subspace 
and the singular loci correspond to singularities of the four-dimensional metric and can be disregarded. 
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where the removed points x A are the poles of the function L. To look at the second 
equation of motion (2.14), it is useful to define L {p) e p (p) such that 

n 

L = ^L {p) . (2.16) 
P =i 

One can then develop (2.14) according to the grading such that 

d*dL m = 
d*dL {2) = 

d*dL (3) = ~[dL w ,[L w ,*dL w ]} 

d*dL w = --[dL w , [L (1 \*dL (2) ]] - -[dL w , [L (2) ,*cLL (1) ]] - -[dL (2) , [L w ,*dL w ]} 
d*dLV = ^[dL {1 \ [L« [L<», [L« *dL«]]]] + l[[L«,dL«], [L« [L w ,*dL w ]]] 

[L (2) ,*dL (2) ]] - ^[dL (2) , [L (1) ,*cLL (2) ]] - \[dn 2 \ [L™ ,*dL w ]] 
o o o 

[LV,*dL™}\ - \{dL m , [L^ML W ]\ - ^[dL (3) , *dL< 1 >]] 
d*dL (6) = ... (2.17) 

and so on and so forth, such that each component L (p) of L can be obtained as the 
solution of a Laplace equation with a source depending of the function L (q) of lower grade 
q < p—1. In other words, one can solve these equations linearly. It should be clear at this 
point that this is a generic property of solutions described by any solvable subalgebra. 
Note that because of the nilpotency (2.1), the system (2.17) stops for a finite degree n. 

The metric and electromagnetic fields defined in four (respectively five) dimensions 
are defined in function of scalars determined as algebraic functions of the component 
functions of L, and 1-forms determined as components of the n valued 1-form W dual to 
the scalar fields 

dW = *VPV~ l = - V ~"~tt ad L fc *dL , (2.18) 

which is well defined on a 

[/(l)xdmi[n] bundle 

over M 3 since the equations of motion 

imply 

d*VPV -1 = 0. (2.19) 

This also provides a Noether charge Q G n associated to any 2-cycle £ of the three- 
dimensional base manifold 

Q E = -j- f +VPV- 1 . (2.20) 

47T l v 
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Regularity. Because we will be interested in composite black hole solutions, we can 
assume that the cycles are characterised by the black holes they are surrounding. In 
particular, since we will study extremal black holes, they will be characterised by poles 
in the function L at point x A located at the horizon of the black holes. By definition of 
the system, the Noether charge will lie in n, and therefore will always be nilpotent. The 
Noether charge associated to a cycle surrounding one single black hole is characterised 
by the pole of *VPV _1 at x = x A . The geometry of the horizon is only modified by the 
existence of the other centres by subleading corrections, and in order for the solution to 
be regular, the associated Noether charge must satisfy at least the constraint that it does 
satisfy for the single-centre solution to be regular. A single-centre extremal solution is 
regular if and only if it can be obtained as the extremal limit of a regular non-extremal 
solution. This requires for instance that Q A is regular in the boundary of the R!j_ x K* 
orbit of a Schwarzschild charge Q Sch = MH [26], 



where (H, E, F) is the SI2 triple inside q which defines the pure gravity truncation. For 
instance, a pure gravity solution would read V = exp(crE) exp(£/H) and the Schwarzschild 
solution would for example be defined such that 



ities requires moreover that the Kaluza-Klein vector oj - the four-dimensional angular 
momentum vector - is a globally defined 1-form on M 3 . The latter will always be defined 
as the specific component of W 

u = ^-TtEW, (2.23) 
k g 

E being still the positive nilpotent element of the pure gravity 5(2 triple inside q. The 
absence of Dirac-Misner string singularities requires therefore 



Q A e <9(R; x K* ■ H) 



dim [iT ■ Q A ] = dim [K* ■ H] , 



(2.21) 




(2.22) 



We define also k 3 such that 2k s = Tr H 2 . The absence of Dirac-Misner string singular 



Tr E Q A — , 



(2.24) 



for all centres. 

On any cycle £(/) surrounding centres x A for A e I, the Noether charge is 




(2.25) 
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A crucial point to be understood, is that although each Q A is required, by regularity, to 
be nilpotent of order 3, The sum of the charges corresponding to different centres, and 
in particular the asymptotic charge of the solution, do not need to be in such a 'regular' 
orbit. Thus, the charge associated to a generic cycle can in general be nilpotent of higher 
degree. Therefore we expect the general solution to admit a total charge 

Q = J2Qa, (2.26) 

A 

which could be any nilpotent element of g, provided the corresponding G-orbit admits an 
intersection with the coset component p. For such a solution, P will be valued in a higher 
order nilpotent orbit on M 3 , although its poles at x A will lie in K* orbits associated to 
regular single-centre solutions. This is depicted in figure 1. As we will see, this remark 
is the key point to obtain interacting non-BPS multi-centre solutions. 

2.2 Parametrizing nilpotent orbits. 

We are therefore interested in classifying the solvable subalgebras n in q in function of 
the nilpotent orbits their generic elements lie in (see [32] for a pedagogical introduction 
to nilpotent orbits). The consistency of the grading with the involution defining K* 
implies that it can always be defined from a particular semi-simple element h (i.e. which 
is diagonalizable in a faithful complex representation) in t* such that for any element 

[h,x {p) ] = 2px (p) . (2.27) 

It appears that the nilpotent orbits are themselves classified in terms of such semi-simple 
generator by mean of the so-called normal triples. Here the factor of two implies that 
the corresponding nilpotent orbit is even, and we will always assume the nilpotent orbits 
to be even. 3 The conjugacy class Gc ■ e with respect to the complex Lie group Gc of 
a nilpotent element e G Qc is entirely determined by the conjugacy class Gc • ft of a 
semi-simple element defining a standard triple (h, e, /) verifying 

[h,e]=2e, [e,f]=h, [h,f] = -2f. (2.28) 

By construction, e is nilpotent, but, from (2.8), we will also require it to lie in the coset 
component p. In this case, h can always be chosen in the t* subalgebra of Q, h G £*, and 

3 This docs not assume any lost of generality because any nilpotent element of an odd nilpotent orbit 
always lies in a solvable algebra associated to an even nilpotent orbit. 
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the triple is then called a normal triple [32]. It will therefore be natural for us to classify 
solvable subalgebras in terms of conjugacy classes of semi-simple elements characterizing 
complex nilpotent orbits under the subgroup A'; C Gc, Kc • e C pc- The classification of 
real orbits requires more work and is of no use for our purpose, so we will not discuss it in 
this paper. Indeed, we would like to straight that also the function of interest L is clearly 
real, the first order system is only characterised by the complex nilpotent orbit it lies in, 
and it will appear in regular solutions that P 'jumps' from one real A"*-orbit from one 
another within the same complex orbit. To summarise, a nilpotent element e G nflp can 
be associated a normal sfc triple (h, e, /), such that its Kc conjugacy class is determined 
by the Kq conjugacy class of the corresponding semi-simple element h G t*. 

Consider the nilpotent element e G p of a normal triple (h, e, /) such that h G t* 
defines an even grading (2.27) of $j. By definition, the positive grade component ^ fc>0 Q (k) 
defines a solvable algebra containing e, and the maximal solvable algebra containing 
Xlfc>oS (fc) * s necessarily a subalgebra of J2k>o8 (k) - Here we will assume that all the 
maximal solvable subalgebras are the strict positive grade component ^2 k>0 Q w associated 
to specific normal triples. This is true for the case of so (4, 4) which we will study in this 
paper, as can be shown by inspection, but we will explain in a forthcoming work that 
this does not exhaust the set of maximal solvable subalgebras for es(8) • 

Now considering a solvable algebra n = ^ fc>0 $j (fe) associated to a normal tripe (h, e, /), 
we would like to parametrize the set of inequivalent embeddings of this algebra inside g. 
The action of K* clearly preserves the graded system of differential equations, and acts 
transitively on the set of embeddings of the normal triple. Physically, the inequivalent 
embeddings of strict normal triples correspond to the inequivalent parametrizations of 
the four- dimensional fields in terms of the set of functions satisfying the graded system of 
differential equations, or in other words, to the inequivalent duality frames the solution 
can be written in. Define KK C K* the maximal compact subgroup of K* and K + (h) C 
K* the parabolic subgroup generated by the Lie algebra Ylk>o ^ {k) associated to h. By 
Iwasawa decomposition, any element u G K* is the product of an element no G KK and 
an element u + G K + (h), u = uq ■ u + . Then, by Sekiguchi lemma [33], a normal triple 
(h, e, /) is conjugate under K* to a strict normal triple (h, e, /) such that h is odd with 
respect to the Cartan involution, i.e. h is Hermitian in a matrix representation. But, 
because uq is invariant with respect to the Cartan involution, any normal triple (h, e, /) 
is conjugate under K + {h) to a strict normal triple (h,e,f). Because K + {h) preserves 
n by definition, the solvable algebra can always be associated to a strict normal triple 
(h,e,f), and so the generator h defining the nilpotent algebra can always be chosen to 
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be odd with respect to the Cartan involution. It follows that the K* orbit of inequivalent 
embeddings of n C g is isomorphically parametrized by the FJT-orbit of the corresponding 
h. Note that for the physical models, KK is always the product of the Ehlers U(l) and the 
maximal compact subgroup K 4 of the four-dimensional duality group G 4 . We have thus 
shown that the inequivalent duality frames corresponding to a solvable system associated 
to a nilpotent orbit are parametrized by the U(l) x K 4 orbit of the generator h defining 
a strict normal triple. In fact we will see that the action of the Ehlers U(l) is physically 
irrelevant, because the inequivalent duality frames it permits to define do not support 
any regular black hole type solutions. The set of 'physically relevant' inequivalent duality 
frames associated to a solvable system of differential equation is therefore parametrized 
by the F^-orbit of an associated strict normal triple generator h. 

In this paper we will restrict ourselves to the example of the STU model, for which 
the three-dimensional scalars parametrize the 'pseudo-quaternionic' symmetric space 

G/K* = SO(A,A)/(SL(2) x% 2 SL{2) x SL{2) x Za SL{2)) , (2.29) 

and so h G ® A =o sl^. Choosing a particular Fa, Fa, Fa basis for each s^ A) , the generator 
h determining a nilpotent orbit is parametrized by four positive half integers b A such that 

h = J2 &A #a • (2.30) 

A 

This parametrization is the key point to obtain the graded system of equations, as we will 
see in the following. Because a nilpotent orbit of K* is always defined as a Lagrangian 
submanifold of a nilpotent orbit of G, it will be relevant to define the conjugacy class 
of h G g by its so-called weighted Dynkin diagram [a^ 2 ]; which evaluates a dominant 
representative of G-h in a Cartan subalgebra on the corresponding simple roots of so(4, 4). 
It is proved that A G {0, 1, 2} (/3 A G {0, 2} for even orbits) [32]. 

The parameters b A associated to a Noether charge Q corresponding to one centre can 
be characterised by the central charges at the horizon as 

6 A (K- 1 gH) = -J_ |Z A .( g ,p)| , (2.31) 
</\h{q,P)\ 

where qA,p A are the electromagnetic charges, and h{q,p) the corresponding quartic in- 
variant, 

h{q,p) = 4Ai<? 2 g 3 - 4<?opW + 4 ^p i+1 p i+2 q i+l q i+2 - (p°q + ^P^i) ■ (2-32) 
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Zo*(q,p) = Z^{q,p) is the central charge at the horizon and Z ijf (q,p) its Kahler derivatives 
in tangent frame. 4 V^QV* describes the pole in P £ p at the centre as a residue, where 
V* is the value of the coset representative at the horizon. Although V diverges at the pole, 
V*~ l QV* is well defined at the horizon because the singular component of L commutes 
with Q by definition of the charge. 

The BPS black holes admit a nilpotent Noether charge such that b A = (2, 0, 0, 0), the 
non-BPS black holes with a vanishing central charge at the horizon correspond to the 
permutations b A = (0,2,0,0), (0,0,2,0), (0,0,0,2) and will not be discussed, 5 and the 
'physical' non-BPS black holes with a non-vanishing central charge at the horizon admit 
a Noether charge such that b A = (1,1, 1,1). 

2.3 Three solvable subalgebras describing all solutions 

There are three solvable subalgebras that contain all the others by permutations of the 
b A, s or as subalgebras, and so we will focus on them. 6 These three algebras have the 
same dimension, and describe three different families of solutions, each of them being 
characterised by 8 harmonic functions. The first is the well-known BPS system, which 
can be interpreted in type IIA as describing D6-D4-D2-D0 bound states. The second is 
the almost-BPS one [9, 10] and describes type IIA bound states where the constituent 
branes of one type have their orientation inverted, as e.g. D6-D4-D2-D0, while the last 
one corresponds to type IIA configurations in which the orientation of half of the branes 
is inverted, as e.g. D6-D4-D2-D0. 

2.3.1 The BPS algebra 

Let us start by the well known example of the BPS Denef system [5, 31]. In that case 
the pertinent subalgebra is determined by b A = (2,0,0,0), i.e. h = 2H . h £ so(4, 4) is 
then characterised by the weighted Dynkin diagram [ /] and defines the following graded 
decomposition of so (4, 4) 

so(4, 4) = l ( - 2) © (2 ® 2 ® 2) ( - 1} © (flli © sl 2 © sl 2 © sl 2 ) (0) © (2 © 2 © 2) (1) © 1 (2) . (2.33) 
The corresponding decomposition of the subalgebra and the coset are respectively 

3 

4 A) = l ( - 2) © (gl, © sl 2 © sl 2 © sl 2 ) (0) © 1 (2) (2.34) 

A=0 

4 i.e. Z t (p,q) = -(t i -?)e-*fl!,(e*2(g J p)). 
They would be BPS in a different truncation of maximal supergravity. 
6 For a discussion of the nilpotent orbits of S'O (4,4) see [34, 20]. 
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and 

3 

so(4,4) e 0s^ A) = 2 ® 2i ® 2 2 ® 2 3 = (2i ® 2 2 ® 2 3 ) ( ~ 1) © (2 X ® 2 2 ® 2 3 ) (1) . (2.35) 

A=0 

The associated solvable subalgebra is the positive grade component 

n BPS (2 X <g> 2 2 © 2 3 ) (1) © 1 (2) , (2.36) 
with generators e^ Q2 G p and e (2) G 5* satisfying 

[ e «ia 2 a3' ^lftiftj = £: ai/3i £ Q2/32 £ a3/33 e( ' ) (2.37) 

all other commutators being trivial. A generic element is nilpotent of order three in all 
three fundamental representations, and it follows trivially from (2.17) that for the Ansatz 

V = exp(L^e« Q2 J, (2.38) 

the equations of motion reduces to the free Laplace equation on L" 1 " 2 " 3 

d*dL aia2a3 = , (2.39) 

and the general solution is determined by 8 harmonic functions on R 3 . So 

JJX\OL-20LZ _ pia2«3 _|_ \ ^ *A ^2 

/ I rr rp \ ^ ' 

A 1 Al 

At each centre, the existence of a regular horizon requires that the associated Noether 
charge 

Q A = U^qT^e^U, , (2.41) 

with 

[ iaia2ai + E e ™ J ' (2 - 42) 

lies in the physical orbit, and each q A must for instance admit a positive quartic invariant 
I4 > ( defined in (2.32) ), that can also be written as 

t j aia 2 a 3 \ _ „aia2a3„/3i/3 2 /33, 7 7i7273 „<5i<5 2 <5 3 ^ n (o ao\ 

J-iK^A I ~ fc ai/3i fc 7i5i fc a2^2 & 72<52 fc a3/33 fc 7353yA Va Va ^ u I^.^OJ 

such that near the horizon x ~ x A 

SL{2) x SL(2) x gL(2) x SL(2) SO (4, 4) 

(SO(2) x SO(2)) x R 2 > (50(2) x S0(2)) x (R8(D © p>) ' 1 ' ] 
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However, it is important to note that the total charge 

Q = E ^ = ^PH^XWa) (E tf^So**) exp(^^eS 7 273 ) (2.45) 

or equivalently g = g A , does not necessarily satisfy such constraint, and one can have 
for example a negative quartic invariant 



tl/3l e 7l5l 5 Q2/32 £ 72 ,5 2 £ Q 3 /33 £ 7 3 ( 53g QlQ2Q3 g /3l/32fe g^^^g <51 ' 5253 < (2.46) 



in which case one has 

gL(2) x SL{2) x gL(2) x gL(2) 50(4,4) 

(50(1, 1) x 50(1, 1)) k R (2) (50(1,1) x 50(1,1)) x (R 8 w ©R( 2 >) l ' J 

in the asymptotic region |x| — > oo, which does not correspond to the asymptotic of any 
regular single centre black hole. The existence of such solution has been exhibited in [35]. 

The dual vector 

dW = -* dL a ^e^ ia2a3 + e ai ^e a2 ^e a ^ :j L a ^ * dl/W»e< 2 > (2.48) 

includes all the magnetic fields and the Kaluza-Klein vector u, although the latter is not 
trivially the component of W along e (2) , but a combination which depends on the explicit 
choice of h G sl^ '. 

The BPS solutions are well-known and have already been treated in the context of 
nilpotent orbits in [30], so we will not come back to them. 

2.3.2 The principal orbit algebra: almost-BPS equations 

The maximal nilpotent orbits of 5O (4,4) are associated to b A = (4,2,2,2) and its 
permutations, which means h = 4Hq + 2Hi + 2H2 + 2H3, i.e. to the weighted Dynkin 
diagram [***], which defines the following graded decomposition of so(4,4) 



so(4,4) = l ( - 5) © l ( - 4) © (3 x l) ( - 3) © (3 x l) ( - 2) © (3 x 1 © l)^© 
(flli © flli © flli © 0(i) (0) © (1 © 3 x 1) (1) © (3 x 1) (2) © (3 x 1) (3) © 1 (4) © 1 (5) . (2.49) 

For a nilpotent element to be really in the maximal orbit and not in a lower one, it should 
correspond to a generic element of the grade 2 component, for which the four elements 
of grade 1 are all non-vanishing. The algebra t* decomposes as 

sl 2 ©sl 2 ©sl 2 ©5t 2 = 1 ( " 4) © (3 x l) ( - 2) © (gl, © © gl, © gi,) m © (3 x 1) (2) © 1 (4) (2.50) 
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and accordingly for the coset component 

2 ©2 1 ©2 2 ©2 3 l ( - 5, ©(3xl) ( - 3) ©(3xl©l) ( ' 1) ©(l©3xl) (1) ©(3xl) (3) ©l (5) . (2.51) 

Such element is then nilpotent of order seven in all the three fundamental representations 
Ri (vector, chiral and antichiral spinor), i.e. Qj = 0, but Q 6 7^ 0. 
The relevant solvable algebra is therefore 

n aSPS = (l©3x © (3 x 1) (2) © (3 x 1) (3) © 1 (4) © 1 (5) . (2.52) 

We will label ef and ef the grade one generators, ef the grade two, e (3) * the grade 
three, e (4) and e (5) the grade four and five. The non- vanishing commutator can be written 
in terms of the symmetric M.+ x invariant tensor Cij k = \£ijk\, as 

[ef , ef ] = ef [ef , ef ] = c m e^ k [e<* ef ] = 8]^ 

[e< 4 \ e^] = e (5) [ef , e (s)J '] = #e (B) . (2.53) 

Only the generators of odd grade are in the coset component, so one can consider the 
Ansatz 

V = exp(-Vef - ITef - ^e (3)i - Me (5) ) . (2.54) 
After some algebra, one obtains the odd 



p = dVe { l) + dlCe™ +(dZi + ^c ilk (K j K k dV - Vk ] dk k )^j e (3)i 

+ (dM + -(VK i dZ' i - 2VZ i dK i + Z^dV) + —c ijk VK j K k (K i dV - VdK*)) e (5) 
V 6 120 J 

(2.55) 

and even component of the Maurer-Cartan form V~ l dV = P + B 

B = hiCdV - VdK^ef + ~ (fcdZi - Z % dk l + ^-c ljk K j K k (PCdV - VdK 1 )) e (4) . 

(2.56) 

The grade one component of the equations of motion (2.14) gives obviously that V 
and K l are harmonic functions: 

d*dV = d*dK i = 0. (2.57) 
Using the harmonicity of V and K l , one reduces the grade three component to 

d*d (Zi + ^c ijk VKiK k j = ^c l3k Vd*d{kiR k ) (2.58) 
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and using all these equations combined, the grade five component gives 




(2.59) 



It follows that with the redefinitions 



a 2 Vfj, 



aZi 



aV 



V + l , K l = K\ 
X -M + Q + l Z t + Vk - l lA VZJC + i 




lv 2 + h (V - l ))c ijk K l K>K k , 



(2.60) 



one obtains the almost-BPS equations, [9, 10], for arbitrary constant a and l\. It should 
be recalled that this system of equations was initially found in a completely different way, 
and it seems quite remarkable that we recover it with our approach. However, one of the 
physical assumptions behind this system being extremality, it is natural that it admits 
a description in terms of nilpotent orbits, which we exhibit here. This provides a group 
theoretical explanation for the graded structure of the almost-BPS system. 

A generic element q°e ( Q ) + q l ef of the grade one component does not commute with 
any generator of © A sl2 A) , and therefore defines a representative of a maximal nilpotent 
orbit. It satisfies = 0. An element with q° = 0, q l ef ) + pie {3)% + p e (5) lies in the BPS 
nilpotent orbit Qi 3 = b A = (2, 0, 0, 0), and can be chosen such that it corresponds to a 
regular single-centre BPS black hole. With one or two vanishing q\ the charge generically 
lies in a nilpotent orbit Qi = and Q,- 4 = for i ^ j. With only one non- vanishing 
q l , the other components can be chosen such that all Qi 3 = 0, but one can then check 
that the corresponding configuration never corresponds to a regular black hole because 
the required reality conditions are not satisfied. 7 For q % = 0, the element is generically 
nilpotent Q 3 = with b A = (1,1,1,1), and can be chosen such that it corresponds to 
a regular single-centre black hole. These requirements can be understood in terms of 
type IIA bound states in a particular duality frame in which q° can be identified to a D6 
charge, q l to D4 charges, pi to D2 charges and po to a DO charge. 

2.3.3 The subregular orbit algebra: composite non-BPS equations 

The next to maximal nilpotent orbits of 5*00(4,4) are associated to b A = (0,2,2,2), 
h = 2Hi + 2H 2 + 2H 3 and its permutations. Its weighted Dynkin diagram is [ 2 »'] and it 

7 By this we mean that the Levi subgroup of the stabilizer of Qi in K* is not compact. 
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defines the following graded decomposition of so(4,4) 

so(4,4) = 2 ( - 3) ©(3xl) ( - 2) ©(3x2) ( - 1) ©( S l 1 © ^ © qI, © sl 2 ) (0) ©(3x2) (1) ffi(3x 1) (2) ©2 (3) . 

(2.61) 

The nilpotent elements correspond to linearly independent doublet of 3-vector of non- 
vanishing cubic invariant, of either equal sign or not. The algebra t* decomposes as 

sl 2 © s\ 2 © sl 2 © sl 2 = (3 x l) ( - 2) © (fa © © © sl 2 ) (0) © (3 x 1) (2) (2.62) 
and accordingly 

2 © 2i © 2 2 © 2 3 = 2 ( ~ 3) © (3 x 2) ( ~ 1) © (3 x 2) (1) © 2 (3) . (2.63) 

Such element satisfies Q { 5 = 0, but ^ 0. 
The relevant solvable subalgebra is 

n nBPS (3 x 2) (1) © (3 x 1) (2) © 2 (3) . (2.64) 

We will label the grade one generators e l a , the grade two generators and the grade 
three generators e a . They satisfy the algebra 

[e«, ej] = ^ k e k [e t , e^j = 5{e a , (2.65) 

with c* jfc = |£ y ' fe |. The relevant Ansatz is then 

V = exp(-K^e l a - M a e a ) . (2.66) 

One computes that 

-P = dK«ei + [dM a - X -c %3k K*E^K^dKl\j e a (2.67) 

and 

B = -\e aP c^ k K^dK^e k (2.68) 
Z J 

and the equations of motion imply that K l a are all harmonic 

d*dK™ = (2.69) 

and 

d * dM a = ^d(e M c ijk K^K^ * dKfj . (2.70) 
17 



This gives us a new system of equations, that we call the composite non-BPS system, 
and we will show in section 5 and 6 that it admits new non-BPS multi-centre interacting 
solutions. 

To finish this section, we note that, as for the almost-BPS case, one can recover 
the physical non-BPS orbit b A = (1, 1, 1, 1) as a subcase: a charge Q = qfe l a + p a e a 
corresponds to a physical non-BPS orbit, verifying = 0, if 

e a mH = ( 2 - 71 ) 

and 

\/i e apP a qP>0. (2.72) 

It is straightforward to check that one can chose several charges Q A satisfying these 
requirement such that Q = Yl Qa defines a generic element of the nilpotent orbit associ- 
ated to b A = (0, 2, 2, 2). The algebra analysis therefore already suggests that interacting 
solutions exist. 

3 From the group theory to the physical system 

In this section we will present all the necessary material required for defining explicit 
four- dimensional solutions from the algebraic solutions we have described in the preced- 
ing section. We will also discuss the expression of the ADM mass in function of the 
asymptotic central charges. 

3.1 Conventions 

Our aim is to describe black hole solutions of the STU model. Since this model has 
already been studied in great details we will only briefly recall some properties of the 
theory, but refer the reader for example to [36, 20, 13] for further details. The STU 
model is a particularly useful truncation of 11-dimensional supergravity, as well as the 
type II supergravity theories. We will sometimes refer to the microstates interpretation 
of black hole solutions in terms of D-branes in type IIA supergravity. The model can be 
understood as Af = 1 supergravity in five dimensions coupled to three vector multiplets 
with the intersection form = \£ijk\- More details about the uplift of the BPS and 
almost-BPS systems can be found in [13], whereas the uplift to five dimensions of the 
new composite non-BPS system will be presented in section 5. 
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We will consider the standard Ansatz for the stationary metric 

ds 2 = -e 2U (dt + u) 2 + e~ w b ilv dx il dx v . (3.1) 

The 8 electromagnetic fields (4 electric A plus 4 magnetic duals A\) of the STU model 
decompose accordingly as 

-2y/2 A Ai b 2 c 3 = Ca 1 b 2 c 3 (dt + uj) + w Ai b 2 c 3 , (3.2) 

where each A^B^Ci runs from 1 to 2 of the corresponding SL(2, R)j duality symmetry 
associated to the complex scalars 



f = di + ie 



1,2,3 



(3.3) 

parametrizing the upper complex half plan U(l)i\SL(2, R)j. As we have already dis- 
cussed, the stationary equations of motion reduce to the ones of Euclidean 3-dimensional 
gravity coupled to a non-linear sigma model defined on G/K* = 50(4, 4)/(50(2, 2) x 
50(2,2)), i.e. (2.13), (2.14): 

Tr P^P U = R ia/ , d * P + [B, *P] = , (3.4) 

where, as we saw in the previous section, P and B are respectively the coset and sub- 
group component of the Maurer-Cartan one-form V~ l dV of the coset representative 
V G 50(4, 4)/(50(2, 2) x 50(2, 2)). The analysis of these equations has been done 
in the previous section. Here, we focus on how the four-dimensional fields (3.1), (3.2), 
(3.3) are encoded inside the three-dimensional coset element V. 

Using the standard convention that capital indices A+, B i: . . . correspond to rigid 
SL(2)i (acting on the left) whereas small ones a,i,bi, ... to local 50(2) (acting on the 
right), we define 50(4,4) as the subgroup of SL(8, R) preserving the metric 

/ £4,R, \ 



\ 









e A 1 B 1 









£ A 2 B2 £ A 3 B3 



(3.5) 



/ 



and the coset representative V G 50(4, 4)/(50(2, 2) x 50(2,2)) as 





( 


A' 

v ai Al 








\ 




V = exp [C AlB2Ca E AlB2 c 3 + <rE] ex P [UB] 







A' 

v ai Al 







(3. 




\ 








{v- x ) B £*{v- 







/ e u v ai Al 


V 



-V„, A 



-U„. Ai 



„— (7„, Di/ 



(^ 1 )B 2 b2 (^ 1 )c 3 C3 
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which must be understood as acting on a vector (A ai , B ai , Cb 2C2 ) ■ Note also that the 
order of the indices for v ai Al is reversed because it is the transverse of v ai Al which is 
involved. A completely explicit form of this conventions is given in the Appendix. 

The v ai Ai are the SO(2)\SL(2, R) representatives 



1 2 

V 2 V 2 




(3.7) 



^a 1 b 2 c 3 define electromagnetic fields, and indices are raised and lowered via S/ Li B i 
e.g. 

/ r~ <~ c. /-D1E2F3 

<UiB 2 c 3 — ^a 1 d 1 £b 2 e 2 £c 3 f 3 (, 
We parametrize the electromagnetic fields in terms of 2 by 4 matrices 



(3i 



(C 



/•111 £112 £121 £122 
£211 £212 £221 £222 




"^3 -^2 -X 1 



-x 



-x L 



(3.9) 



and in terms of 4 by 2 matrices from (3i 



(a 



/ Cm C211 \ 

C112 C212 

C121 C221 

V C122 C222 / 



f-x° x 1 \ 

X 3 ~tp2 

X 2 -fa 

\ -1p! -Tpo ) 



(3.10) 



Triality is realised as the permutation of the three indices % = 1, 2, 3 in this basis. 
The involution defining the K* subgroup is defined from the twisted transpose 



/ as 

5al 

\ 



\ 





/ 4 1 



















\ 



(3.11) 



~0cl ) 



such that V^" = V 1 iff V G K* . The relation to the representative exp(— L) in the 
symmetric gauge (2.7) discussed in the previous section is 



exp(-2L) = VV X 



x e~ 2U M A ^a + le- 2U M c ^( A ^ E ^ ClD2E3 



x 
x 



e ~ 2U M AlBl 

e M \DiB2C3 



\ 



(3.12) 
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where 



M 



AiBi 



-2<t>i 



+ e 2(pi a 2 e 2 ^ai 



(3.13) 



and we only wrote the elements of the middle column because they are enough to de- 
termine all the fields, when using triality. We recall that the interest of looking at W* 
instead of V is that VV* is, by construction, invariant under [SX(2,IR)] 4 gauge transfor- 
mations. From (3.12) it is then easy to decode the three-dimensional matrix to obtain 
the four-dimensional scalars. To obtain the full four- dimensional fields, one finally needs 
to dualize some of them to obtain the expressions of the vectors. 

From the equations of motion, the Kaluza-Klein vector u is defined by 

1 



du = Tr E dW = *e 



-4U 



da _ -C A ^dC Al B 2 C 3 



(3.14) 



and the magnetic vectors by 
1 



dw 



Tr E AlB2C3 dW 



e- 2U M^ Dl M£ E2 M£ F2 * dC D ^ -e~ w * (da-^^dCn^UrB^ 



e- 2U M^ D M^ E2 M^ F2 *d( D ^ - ( AlB2 c 3 du 



(3.15) 



Note that the electromagnetic field strength are then manifestly twisted self-dual 

-2V2F AlB2 c 3 = dU lB2 c 3A {dt + u)+e- 2U Mxl Dl M^ E2 Mc 2 \^dC D ^ . (3.16) 

This ends the parametrization of the physical fields in terms of the three-dimensional 
coset element V. 

In order to obtain solutions associated to a given solvable system as discussed in 
the previous sections, we must now provide an explicit basis of function L lying in cor- 
responding solvable algebra n D p. In this aim, we define a particular basis of Cartan 
generators of the subalgebra Q) A sl^ (according to the matrix notation) as 



where Ua are the four specific 4 by 4 matrices 




(3.17) 



U A 








o 



°A 
1 _ 
2 



2 °A 



/> 3 I - r) 2 
°A 2 °A 








(3.18) 
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The generic elements X of p are defined similarly as (3.22) below. One can easily solve 

[H A ,Xp]=p A X fi , (3.19) 

for /?a = (±1, ±1, ±1, ±1), which defines a basis for the 16 elements of p. The explicit 
form of the elements Xp is given in the appendix. The functions L of (2.7) in the positive 
grade components of p associated to a solvable system are then simply obtained in this 
basis as 

• l = ^2 lPx p in the BPS s y stem ; 

/3o=l 

• L = L^Xp in the almost-BPS system; 

2/3o+Ei ft > 1 

• L = L^Xp in the composite non-BPS system; 

£ift>i 

respectively. 

As explained in section 2.1, the most general solutions of these types are generated by 
the action of KK = [?7(1)] 4 . The general Cartan basis H\ can be obtained by conjugating 
(3.18) by a general [£7(1)] 4 rotation, given explicitly in the appendix, such that H\(a\) 
is a function of the phase «a only, by property of the ® A s[2 A) algebra. Nevertheless 
it will be easier to first compute the explicit basis Xp for specific H\ and then rotate 
the associated function L = ^2 L^Xp with respect to KK, rather than to compute the 
explicit basis for the general H^(a\) from scratch. 



3.2 Charges and central charges 

We are now going to describe how the central charges of the solution can easily be 
extracted from the coset representative V, and how the mass formula can be determined 
from the algebraic conditions satisfied by the solutions. This will permit to generalise 
the standard BPS formula M — \Z\ to the non-BPS systems. 

Considering any space-like cycle £ embedded in M 4 through l, orthogonal to the 
time-like vector d t , the associated electromagnetic charge is 

<?AiB 2 C 3 |£ = 2^J^ f^MB^Cs =_ 47rv /2 J d{ w AiB 2 C 3 + ^CaiB 2 C 3 ) = ^/| Tr ^A 1 B 2 C a Q\^ 

(3.20) 
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where we used the absence of Dirac-Misner string 



Tr E Q\y, = 



(3.21) 



to show that the contribution from CaiB 2 c 3 w drops out. 

It is straightforward to compute that the total Noether charge Q can be rotated to 
the coset component using the asymptotic value of the scalar fields V 



\ 



NS a d l 

Mi 



Z' 



a\e 2 h 



\ 



~ z a\e 2 fz 



z d\b 2 c 



(3.22) 



°b 2 °c 3 % a c 3 / 



defined as a matrix acting on a vector (A dl , B dl , C e2 f 2 ). Here, the (2,3) component 
(— z aie2 f 3 ) and (3,1) (— z dlb2C3 ) can be understood as matrices, as minus the transpose of 
(z 0162 ^ 3 ) and (zd^c?) defined as in (3. 9), (3. 10). M is the ADM mass of space-time, the 
total NUT charge N is assumed to vanish, a a . bi are symmetric traceless matrices 



(3.23) 



associated to the momenta of the scalar fields in the asymptotic region 

r 2 d r f 



II 



lim 

r— >oo f l — t l 



and 



aib 2 c 3 _ ai b 2 c 3 A 1 B 2 C 3 



(3.24) 



(3.25) 



defines the [SX(2,IR)] 3 covariant combination of the asymptotic central charge and its 
asymptotic Kahler derivatives, which we will call the asymptotic 'central charges'. The 
asymptotic 'central charges' are defined in terms of z aib2C3 as 



Z 



-l ) 0,1 + a2 + a 3 Z a i a 2 a 3 



(j\ a i ( _j\ a i+i+ a i+2 z aia2a 3 



at ,02,03 



We have 



Qo 



(3.26) 



(3.27) 



and the Z^'s are obtained by replacing by its complex conjugate in the holomorphic 
component. 
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It will be illuminating to discuss the action of KK = [U(l)} 4 on the asymptotic 
momenta we just defined. These rotations act on V ~ 1 QV given in (3.22) such that 

M — iN = ei (Q0+cn+a2+a3 )(M — iN) Z = e i(- Q o+ai+«2+a3) % 

~ 4 (3-28) 

JJ. = g2 ( — a °~ a i+ a i+l+ a j+2)TJ. 2j = gi( _Q 0+ai-Oi+l— ai+2)g. 

Considering a solution associated to a given nilpotent orbit, V ~ 1 QV lies by definition in 
the closure of this nilpotent orbit. It follows that it lies in the positive grade component 
defined by a generator h = b A H\(ctA) for some phases a\. We now have to specify the 
b A i. e. the orbit we are in, to obtain an explicit relation. One can already anticipate that 
the number of phases it depends on is simply the number of nonzero b A labelling the 
orbit. 

For the BPS case, it is well known that one has 8 

M-iN = e~ ia °Z, U i = e ia °Z i . (3.29) 

This translates the fact that only b° ^ in the BPS orbit. Requiring the absence of 
NUT charge determines the phase ao as the one of the asymptotic central charge such 
that M = \Z\. 

We now turn to the almost-BPS orbit: in order to obtain the general form of the 
asymptotic momenta in function of the asymptotic central charges for the non-BPS sys- 
tems, it is convenient to first compute the constraint at a a = 0, and then obtain the 
general solution by applying the transformations (3.28). The property that the almost- 
BPS system can be obtained from the BPS system by replacing a D6 charge by D6, 
illustrates into the fact that the ADM mass and asymptotic scalar momenta are given by 



M-iN = Z - -Re 
2 



J2 Zi] , Ui = Zi- ^Re [Z + J2 Z \ ■ ( 3 - 3 °) 



.i 



With respect to the BPS case, the new term comes from the presence of the D6 instead 
of the D6. One can indeed check, using (A. 5), that for a pure D6 charge, one obtains 
M = —Z. Reintroducing the angles, we see that in this case the ADM mass is not 
determined by the asymptotic 'central charges' only, but depends on three phases ccj 
such that 



8 Here we shift all phases by ^ with respect to the convention (3.18) to avoid the presence of extra 
i factors. 
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where «o is determined in function of the a^s such that this expression is real positive, 
and in particular such that the NUT charge vanishes. 

The composite non-BPS system asymptotic momenta can themselves be obtained 
starting from the almost-BPS system by substituting D4 charges to the D4 charges as 



M-iN 



1 



Z - -R 



[ z + E^HE Im [ z -z< 



(3.32) 



and 



Zi Re 

2 



Z 



-Z + Zi 



H — Im 

2 



E^ 



z 3 +y^. z k 



(3.33) 



It is remarkable that although the expression of the ADM mass and asymptotic scalar 
momenta are rather complicated deformations of the BPS expressions, they end up being 
pretty simple. Restoring the dependency in the phases one gets 



M-. 



n, 



l Z + e 



i(ai—a i+ i—a i+ 2 



H+l 



-iai+3 



z, 



+2 



(3.34) 
(3.35) 



where the phases ccj are again assumed to be chosen such that the ADM mass is real 
positive. Note that the dependence in ao drops out, which comes from the fact that 
b° = in the composite non-BPS case, and thus the associated rotation leaves invariant 
h = 2 Yli Hi{oti). It is remarkable that the ADM mass reduces to the ADM mass formula 
computed in [17] within a specific Ansatz up to a four- dimensional duality transformation 
parametrized by the three phases CKj. It follows that the ADM mass of a bound state of 
non-BPS black holes is defined by the generalised 'fake superpotential' defined in [21], for 
specific values of the 'auxiliary fields' associated to the flat directions which will depend 
in general on the specific interior structure of the solution. Note that the 'auxiliary fields' 
are extremum of the generalised 'fake superpotential' if and only if there is no interactions 
between the centres. 

To finish, let us discuss the case of single-centre non-BPS black holes, that are par- 
ticular solutions of both the almost-BPS and the composite non-BPS system. They 
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furthermore satisfy that 

I m [ e |(-«0+ai-ai+l-a i+2 ) Z .j = j m ^(- ao+ a 1+ a 2 +a s ) ^ ^ (3.36) 

which corresponds to the absence of D4 or D4 charges respectively. All the phases oa's 
are then determined in function of the asymptotic 'central charges' as in the BPS case, in 
accordance with the no-hair theorem. One computes straightforwardly that the duality 
rotation 

Z = e iE ! (« ! -iE A «A+^) Z ^ Z . = e i(a 4 -i Ea«a+^) z . ^ (3.37) 

provides a solution to the non-standard diagonalization problem which defines the 'fake 
superpotential' [20] (with a = \ XIa^a 1X1 (2-68) of [20]), such that (3.34) indeed reduces 
to the asymptotic value of the 'fake superpotential'. The equation of the asymptotic 
scalar momenta implies the existence of two flat directions 

Im[ e i(- Q °- Q '+ a *+i+ a *+2)nJ = -lm[ei { ~ ao+ai+a2+a:i) Z] . (3.38) 

It is important to understand that the action of KK = [£7(1)] 4 on the solvable systems 
themselves is very similar to its action on the asymptotic momenta. From the stabilizer of 
the defining semi-simple element h = b A H\ one obtains similarly that the BPS system is 
parametrized by one single phase a , the almost-BPS system by the four phases «a, and 
the composite non-BPS system by the three phases The discussion of the last para- 
graph suggests that there is always one combination of the parametrizing phases which 
does not permit to obtain new interesting systems, as far as asymptotically Minkowski 
composite black hole solutions are concerned. Indeed, the BPS system is known to be 
unique. We will see in the following that the Ehlers rotation of the almost-BPS system can 
also be reabsorbed in a reparametrization of the solution, excepted for a particular value 
of the phase which gives rise to a system which does not admit regular composite black 
hole solutions. It follows that the general almost-BPS system is in fact parametrized by 
three phases associated to the compact subgroup of the four-dimensional duality group. 
The Ehlers U(l) acts already trivially on the composite non-BPS system, and there is no 
further restriction in this case, such that the three phases then parametrize inequivalent 
classes of solutions. 

4 Almost BPS solutions 

We have seen in section 2.3.2 that the equations associated to the maximal orbit of 
so(4,4) are exactly the almost-BPS equations [9, 10]. Here, using the tools from the 
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previous section, we will relate the system of equations, obtained from a purely algebraic 
point of view, to physical solutions of the STU model. We will first show how to recover 
all the known almost-BPS solutions [10, 11, 12], including their "generalised version" 
[13, 14], and show that the Ehlers rotation does not produce any new solutions than the 
one obtained by the four-dimensional dualities. 

4.1 Recovering the almost-BPS usual solutions 

We now want to decode the four-dimensional fields out of M. = Vf]V t i]. This is of 
course not a one-to-one correspondence, because the fields from the three-dimensional 
perspective are duality invariants. One must thus choose a duality frame in order to write 
an explicit solution. It is interesting to choose the one where V correspond to a D6-charge 
and K % to D4-charges. Indeed, this is the duality frame where the first almost-BPS 
solutions have been found, and is therefore well-known [10]. This corresponds naturally 
to the parametrization given in the previous section and in the appendix (3.17), (A. 5). 
The general solution should include all the dependences in the KK angles, but, for clarity, 
we will not consider it here. It should however be clear that the three free phases of KK 
are sufficient to obtain the whole class of almost-BPS solutions discovered in [13], and 
this in a straightforward way. 

In order to recover the almost-BPS solution in its original parametrization, it is 
interesting to use the freedom of redefinition of (2.60). We here assume for convenience 
the constants in (2.60) to be a = —2 and lo = 1% — —2: 



With these choices, After having chosen this particular duality frame, the four- 




(4.1) 
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dimensional scalar fields can be extracted from V using (3.12): 



e 



e - 4U = VZ 1 Z 2 Z 3 - V 2 fj? 

2<f>i 



Zi 



X° = e 4 ^V, (4-2) 
x * = l + e 4U V(-Z l+1 Z l+2 + VlCfj) , 

= -1 + e 4U ( Z 1 Z 2 Z 3 - V> ^ Z l K i + Z l Z 3 K l K 3 J , 

^ = e iU V ({Z i + VK i+1 K i+2 )n- Zi^tfzA , 



a=—^-^(-2 + V)fji + f i^2(Zi + VK i+l K i+2 ) - ]T ICZ^ 

i i^j 



This describes exactly the solutions of [10, 11] 



9 



Vectors. For completeness, one also have to check that the vectors are the same. These 
ones are given by the equations (3.14), (3.15). Having chosen a particular duality frame, 
it will be more convenient to write the electro magnetic fields as two four-vectors. The 
map from the 2 x 2 x 2 to two four- vectors is given by (3.9), with the equivalent for the 
vectors w A , v a- 10 

A A = X A {dt + u) + w A , (4.3) 
A A = tp A (dt + u) +v A , 

The to equation rewrites 

★ dco = e 4U (da - l - ( X A #a - i>Kd X k ) ) . (4.4) 

Plugging in the values for our system (4.2), this equation greatly simplifies, to give 

*du) = dtyn)-VZ i dK i . (4.5) 

9 Because these fields do not strictly speaking define tensors, we have not used Einstein summation 
convention in these equations. We recall that these solutions have been first written in a five-dimensional 
language. The map to its four-dimensional version is done in [13]. 

10 Note that we do not introduce the — 2\/2 factor in the definition of the four-dimensional vectors as 
in (3.2) in order to facilitate the comparison with [10, 11]. 
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One can also see that using (4.2), the equations for the electromagnetic vectors (3.15) 
simplify drastically: 



★ dw° = dV , *dv = ZidK* - K l dZ { + {Vd(K l K 2 K z ) - K l K 2 K*dV) , (4.6) 
M = - ★ dw - VdK { + ICdV , *dvi = dZi + ^- (K j K k dV - Vd(K j K k )) . 

Up to a small gauge choice, this is exactly the equations for the vectors of [13] . Therefore 
one completely recovers the known almost-BPS solutions. 

4.2 Ehlers rotation 

In [13], the authors have used the four- dimensional duality group SX(2,R) 3 of the STU 
model to produce new solutions. As we have already seen, the transformation that can 
be done with this group are part of KK. However, in three dimensions, there is another 
£7(1), coming from the Ehlers reduction along time. In our framework, it is thus clear 
that one can reobtain all solutions of [13], by acting with KK on VV*. In principle one 
could also obtain new solutions with the extra Ehlers U(l) action, but we will see that 
this is not the advocated in the previous section. 

To show this, we perform a rotation with angles 11 «o = oil = 0:2 = «3 = a/2 

(yyty = R( a ) VV % R{a) T . (4.7) 

The expressions for the physical fields are then a priori a complicated mixture of all 
involved functions, and a factors. However, one can do the following redefinitions 

V=-(v + l ), K i = -k i + mna, 
a \ / c 

Zi + hv + 3l )K i+1 K i+2 ) +k + V 2^7 tan a, (4.8) 



a c 2 



\ i i 

-2 tan a (l + ^ K { K j - tan a ^ % - tan a V + 2 tan 2 a ^ K { + m + 

i<j i i 

These redefinitions seem to be rather complicated, but the point is that they are just a 
generalisation of (2.60), where we had an extra harmonic function in the Z^s and in M, 



n The additional rotations along the SL(2, R)j, i = 1,2,3 are done by convenience, and do not change 
the argumentation. The precise parametrization of the rotation is given in the appendix. 
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and an additive and a multiplicative factor to the Zfj's, but one can easily check that 
this is still compatible with the almost-BPS system of equations (2. 58), (2. 59), with the 
constants being given by 

2 r r 2 . sina 

a = — , c = cos a , l = k = , m = 4 — — . (4.9) 

cos^a cos a cos z a 

These redefinitions degenerate for a = tt/2, but one can check in this case that there are 

no physical solutions. Then e~ 4U and the scalar fields are exactly given by (4.2) : 

e~ 4U = VZ 1 Z 2 Z 3 -V 2 ^ 2 , e -^ = e —, * = (4.10) 

The electromagnetic fields are also, up to a constant, the same as in (4.2) 

X° = -sina + e 4C V 2 /z, (4.11) 
x * = C os a + e 4U V(-Z i+1 Z l+2 + VK'/i) , 



= _ l + srn^ + e4U f zzz z^ + yV Z,Z,K : lA , 

cos a \ I 

i i<j 

ifji = sin a + e 4U V I {Z t + VK i+1 K l+2 )fi - K j Z 3 J . 

So we are in the same class of solutions as before. Note that, because of the different 
constants in the electromagnetic fields, the explicit expression of a differs from (4.2), but 
they coincide up to a gauge transformation, and are therefore physically equivalent. We 
conclude that the Ehlers rotation does not produce any new solutions, but led us in the 
class already spanned by the four- dimensional duality group. 



5 Composite non-BPS solutions 

We now turn to the study of the composite non-BPS class of solutions unveiled in section 
2.3.3. Because this class was not known so far, we present it in more details. We will 
focus on the four- dimensional solution, although the uplift to five dimensions will also be 
shortly discussed. In the next section, we will exhibit an explicit two-centre solution in 
this class, and study its physical properties. 



5.1 Equations of motion 

This system represents solutions where half of the branes have been inverted. The major 
difference with the almost-BPS system is that our system is now "symmetric" in terms 
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of the two triplets of functions Kt a=1 and Ki a=2 , whereas the latter was not in terms of 
K % and Zj. It is however convenient to break this manifest symmetry by redefining the 
following functions 

Kr 1 = 2(L i + l), 

Kt =2 = K^ (5.1) 
M a=1 = 8M+Uj2 K ' + J2 K ^- 3 J2 K i L i+i L i+* ) > 

M a = 2 = 2(V + 1) + 1 1 2 K t K 3 - 3 K i+1 K l+2 L t J . 

With these redefinitions, the system of equations becomes 
d -k d Ki = , 

d*dL i = 0, (5.2) 
d*dM = -& k d(LiLj * dftTjfe) , 

d*dV= -c ijk d(Li * d(KjK k ) - KjK k * dL<) . 

5.2 Reconstructing the solution 
5.2.1 Four dimensional point of view 

As we have done for the almost-BPS solutions in the previous section, we want here 
to reconstruct the four- dimensional solution from the matrix. We will not present 
the most general solution, but restrict, for the sake of clarity, to one particular duality 
frame. In order to compare easily with the results of the almost-BPS class of solutions, 
we uplift the solution in the duality frame where (V, Ki, Li, M) correspond respectively 
to the (D6,D4,D2,D0)-charges. In particular, for Ki = 0, one falls back into the physical 
non-BPS orbit, describing single-centre D6-D2-D2-D2 black holes. In our duality frame, 
the redefinitions (5.1) are such that e~ iU takes the simple usual expression 

e ~ iU = VL X L 2 L Z - M 2 . (5.3) 

Introducing the quantities 

Ti = VLi - 2KiM + K 2 L l+1 L l+2 , (5.4) 
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the scalar fields are given by 

t i = KjLj+iLj+2 -M + ie~ 2U ^ ^ 

The electromagnetic fields and a are finally given by 

x ° = e AU (-MV + (VL^L^Ki - MK l+1 K i+2 U) + L X L 2 UK X K 2 K^ , 

i 

x * = -1 + e 4C/ (M(K i+1 L i+2 + K i+2 L i+1 ) - L l+1 L l+2 (V + UK l+l K i+2 )) , 

ij = -l-e 4U L 1 L 2 L 3 , (5.6) 

^ = - e 4C/ L,(-M + KiL i+1 L i+2 ) , 

a = — ( K -M(2 + V + Y^)+Yl (VLi+iU+iKi - MK t+1 K l+2 L t ) 

i i 

i 

Vectors. Having chosen a duality frame, one can compute the explicit equations for 
the vectors. The general equations are given in (3. 14), (3. 15). Using the explicit values 
for our system (5.3)- (5.6), these equations greatly simplify, to give 

* du = dM - ^c ljk LiLjdK k (5.7) 



for u), and 



* dw° = -dV + ^(LidiKjKk) - KjKkdLi) , 
i ,dw i = *du + c ijk (KjdL k - L k dKj) , (5.1 

*dv = -kdu = dM — ^c ljk LiLjdK k , 
-kdvj = dL 



i ■ 



for the electromagnetic vectors. One can check, as it should be, that the compatibility 
equations are nothing but the scalar equations (5.2). 

5.2.2 Uplift to five dimensions 

It is interesting to uplift the solution to five dimensions. In this framework, the fields 
are simply the five dimensional metric ds^, the three electromagnetic one- forms A\ and 
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three real scalars X 1 , % = 1,2,3 verifying X 1 X 2 X 3 = 1. They are related to the four- 
dimensional fields through 



d S 2 = e 4*/3 (# _ A ) 2 + e -2*/3 ds 2 > 

A\ = A 1 + Re(f ) (dip -A ), (5.9) 

X i = e -2*/3l m (fi) , 



where ip is the fifth coordinate, and we defined $ by 



e 



2 * - Im^rin^Im^ 3 ) . (5.10) 



Applied explicitly to the solution of the last subsection, this gives the five-dimensional 
metric 

ds\ = -!L(dt + u- e" 4l/ |W " ™°)) + - ^°) 2 + Tdsl , (5.11) 

where T is defined as T = (TiT 2 T 3 ) 1//3 , and H by 

H={y- C —K,K 3 L k ^ + AK^K^M - C —K i L j L^j . (5.12) 
The one-forms are given by 

A\ = dt--(v+ (LiK i+ iK i+2 - KiK i+1 L i+2 - KiK l+2 L l+1 )^ (dt + u) 

i KiL i+ iL i+ 2 — M , x , . 

-I- w — cu H — [dip — w ) . (5.13) 

J- i 

In this duality frame, the presence of the dip — w° term in the metric implies a non-trivial 
fibration of the associated circle over the four-dimensional base in the presence of a non- 
vanishing D6 charge. In particular, fixing the D6 charge to 1 and choosing appropriate 
boundary conditions such that T 2 and H both scale as r~ 2 in the asymptotic region 
should permit to define asymptoticaly Minkowki solutions. It should be interesting to 
investigate this question. 



5.3 Rotation along time 

We show here that, as expected from the fact that b° = for the composite non-BPS 
system, the Ehlers rotation does not produce any new solutions, and thus that the whole 
class of composite non-BPS solutions is parametrized by the four- dimensional duality 
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group. In fact, the rotation along the time direction uses the symmetry in a = 1,2. To 
show this, starting from the general given in (3.12), we perform a rotation, with 
general parameter a/2 along the U(1)q in KK, together with a rotation along each of the 
three U(l)i with parameter —a/2 

(yyt)' = R( a )VV x R(a) T . (5.14) 

Then, if one redefines 

K\ = 2((1 + cos a + K { sin a) , 
K i = 2 [Ki cos a — (1 + L») sin a) , 



+ ^4(1 + V) + | ^ - 3 £ LiK^K^ j sin 
M 2 = ^2(1 + V) + ~ KiKj - 3 £ Li^i+i^J j cos a 
- |^4M + ~ ^ ^ + KiLj - 6 £ L l+1 L l+2 K^ j 



cos a (5.15) 



sm a 



one first can show that the new functions still solve the system of equation (5.2), and 
that one recovers exactly the fields given in the previous subsection. This proves that the 
rotation along the Ehlers SL(2, M) just transforms a solution into another of the same 
class. 



6 Solving the composite non-BPS system 

In the previous section, we have written the four-, and five-dimensional Ansatze associ- 
ated to the new system of equations (5.2), and we will now solve the equations explicitly. 
As we have explained previously, for a solution to be into this new orbit, and be regu- 
lar, it has to have more than one centre. Indeed, the single-centre solutions are forced 
by regularity to be in the BPS or the physical non-BPS orbit. We will not discuss the 
single-centre non-BPS solution in this paper, and refer to [19, 10] for an explicit exposi- 
tion of the latter. For the sake of simplicity, we will only construct here an axisymmetric 
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two-centre solution. 12 Similarly as for the almost-BPS class of solutions, we are going to 
see that the system admits regular solutions that do not belong to any previously known 
subsystem. Although regularity implies that the coset momenta P falls back into the 
physical non-BPS orbit at each centre, it belongs to the higher order subregular nilpotent 
orbit at a generic point of the three-dimensional base subspace. The centres are then in- 
teracting, in the sense that their electromagnetic charges produce an angular momentum 
and the distance between the centres is fixed. 



6.1 The solution 



One parametrizes the three dimensional base space in spherical coordinates (r, 8, <p) and 
take the two centres to be at r = for the first one and along the positive z axis (9 = 0) 
at a distance R from the origin for the second one. We denote the polar coordinates 
centred at the second centre position as (£,#£). Their relation to the polar coordinates 
(r, 9) centred at the origin is: 

r cos 9 — R 



£ = Vr 2 + R 2 -2rRcos9 



cos 9s 



(6.1) 



Looking at a general two-centre solution, we take the Lj's and K^s to be of the form 

(6.2) 



r 2j r L 



We now have to solve the equations (5.2) for V and M. Let's first look at the equation 
for M. From (6.2), it can be rewritten as 



c ijk 

d*dM = d 

2 



A *d 


dk dk 






r S 



r S r 2 S 2 rS 

It will be convenient to decompose these linear equations into the several pieces 



. (6.3) 



d*dM 1 = d(-) A ★<£(-) , d*dM 2 = d(— ) A*d(-) , d*dM 3 = d[-) A *d(—) , 
d*dM A = d(\) A*d(^-) , d * dM 5 = d(—) A *d(~) (6.4) 

and symmetric expressions in r -h- S. By direct computation, the solutions to these 
equations are 



2 3r 3 



1 



2rS 



M A 



cos 9 



Ms 



COS 

2i?r 2 ' 



(6.5) 



12 Note that a two-centre non-BPS solution is not axisymmetric if the intrinsic angular mementa of 
the black holes are not parallel to the line which joins them. 
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We obtain symmetric expressions in r and E by exchanging (r, R, cos 9) with 
(E, — -R, cos#s)- Finally, including the allowed harmonic contribution for M, the solution 
is 



m m cos 9 _cos#£ 
M = m H h — + a — — + a——— 



Jjk 



UQjdk hQjdk QiQjdk QiQjdk h(Qjdk + Qjdk) 



3r 3 

+QiQjdk~j^=; ~ QiQjdk R ^ - 

Because one can rewrite the equation for V as 

c ijh 



3E 3 



rE 

+ QiQjd, 



(6.6) 



COS I 



i?E 2 



d*d(V - —KiKjLk) = -c iik d{KiKj) A*dL h 



(6.7) 



the structure is exactly the same as for M, multiplying the terms by —2 and exchanging 
Ki and Lj. Hence, the solution for V is 

Jjk 

T 

Jjk 



V=—L i K j K k + h H h — + /? — g- + /J 2 

KdjQk . kdjQk . didjQk didjQk k i (d j Q k + djQ k ) 



r* is 3r 3 3E 3 rE 

+didjQk - — — didjQk — didjQk 2 — h didjQk- 



3E 3 



(6.8) 



/.'/•>: ' /.v>: ' /,'r 2 ' ' ' re 2 

To obtain the full solution, we finally need to solve (5.7) and (5.8) for the vectors u, 
w A and v\. Clearly, the terms Mi and M 2 in M involving only one centre will give no 
contribution to u. The different contribution in the equation are then 



-k duii = d( — — 
V2rE 



r h 



•kdui2 = d 
cos 9s 



9\ 



RrZ. 



1 

r z E 



(6.9) 



•kduiz = d 

Again by a brute force resolution, the solutions are 
R cos 9 — r 



2Rr 2 



1 A 1 
—d~. 

rl-i r 



U)i 



2i?E 



-d(j> , C02 



sin 2 9 sin 2 



i?E 



This leads to the full solution for ui 



oj ■ 



Jjk 



m 



liljdk ) cos 9 + [m 



Jjk 



Uljdk ) cos^e — a 



2i?E r 
sin 2 # 



(6.10) 



.sin 2 0, 



a- 



iifc/i /p; , „ -r ,r - _Rcos6> ^^~sin 2 # 

+c tJ *(Z i (Q J -d* - Qi4) — ^ QiQjdkT^ 

r sin 2 6* 



2i?E 

~ ~ r sin 2 ~ sin 2 9 
+QiQjdk QiQi k 



(6.11) 



2i?E 



2i?E 5 
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The vectors v\ are then easy to find. Indeed, from (5.8) vq is just equal to u, and the 

(6.12) 



Vi's are the harmonic duals of the Lj's functions. Then we have 



V = UJ, 

Vi = Qi cos 9 + Qi cos 0£ , % — 1, 2, 3 . 

As for V, finding the solution for the u> A does not involve any new terms, and the solutions 
are given by 

,sin 2 # —sin 2 6*5 
r 

r — R cos 9 , , ~ sin 2 9 



w 



(Qs + c ijk kikjQ k ) cos 9- (Q 6 + c ijk kikjQ k ) cos # s + /?^— + /?- 



+2& k \k i (d j Q k -d j Q k 
r sin 2 



-didjQk- 



2i?S 
sin 2 # 



didjQk 



RT? 



- didjQk—— + c/i^Qfc 



RE 
r sin 2 

2RT? 



(6.13) 



u> =u — c 



(kjQk - ljd k ) cos9 + (kjQ k - ljd k ) cos9- E + (d j Q k - djQ h 



r — R cos 9 



6.2 Regularity and physical properties 

The solution written in the previous subsection solves the equations of motion, but does 
not a priori satisfies the regularity conditions, necessary for the solution to be physical. 
We investigate this question in this subsection. One notes that the regularity conditions 
are exactly the same in four and five dimensions. However, the solution we construct 
being asymptotically Minkowski in four dimensions, the physical charges are the four- 
dimensional ones. 



Dirac— Misner strings. We will first implement the conditions for the absence of 
Dirac-Misner strings singularities. This requires u to vanish for sin 9 = 0. Because 
of the presence of the two centres, this imposes three conditions, 



Ue=iT — , W0=O,r<R — , &6=0,r 



>R 







They can be solved by affecting the values of k, m and m to 



k = [Qjdk — Qjdk 
Solving for k, one can rewrite it as 



m 



c ijk 

-k + —hljd k 
2 



m = k + —hljdk 



m + m 



J,jk 



—m + 



-Ijljdk 



c ijk 

—' 

c ijk 

2R l 
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-hlj(d k + 4) 



h ( Qjdk — Qjdk 



(6.14) 

(6.15) 

(6.16) 
(6.17) 



The second equation is similar to the non-BPS version of the integrability equation, or 
bubble equation [5, 31, 37, 38, 39]. However, it should not be understood as such. Indeed, 
using (5.8) one sees that this condition is equivalent as to the vanishing of the DO charge 
of the two black holes. Nevertheless, we will show in the following that the distance 
between the two centres is determined in function of the electromagnetic charges of the 
two black hole constituents and the asymptotic value of the scalar fields. The property 
that this constraint does not directly follow from the absence of Dirac-Misner string 
singularities is peculiar to the simple duality frame we are considering, and we expect it 
would be fixed in this way as in the BPS case if we were considering the most general 
solution, function of the three phases ctj defining the duality frame. 

Horizon regularity. We then have to make sure that each of the centre is a regular 
black hole. We perform this analysis at r — > 0, the result at £ — > being exactly 
equivalent. The horizon area is given by 



A = J d9d^^r 2 (e- AU r 2 sm 2 e-uj 2 ) (6.18) 



and it is therefore necessary for this quantity to be finite as r — > 0. One has 

r^o sin 2 6 

u = —a \-U(r), (6.19) 

r 

which implies that the contribution from u 2 is always finite and well-behaved. On the 
contrary, as r — > 0, e~ AU diverges in general as 1/r 6 : 



-4U r->0 



^ {^QxQ2Q^ ik Qidjd k ) - ^ h QiQjd k ) 2 J + 0(r~ 5 ) . (6.20) 



This term vanishes if 



— — — — — — 7 (621) 

Qi Q2 Q3 

for some arbitrary constant 7. With this condition, e~ 4U still diverges as 1/r 5 : 

e -4U r^O Q1Q2Q3W - 27a) COsfl + Q ^_ 4 j ^ 

which imposes us to take 

/3 = 2 7 a. (6.23) 

This conditions are now enough to ensure that the black hole located at r = has a 
finite entropy, and one shows that they are necessary. The argumentation is exactly the 
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same for the centre at £ = 0, and therefore we also imposes 

di d 2 d 3 „ 

— = — = — = 7, (6.24) 

Qi Q2 Q3 

/3 = 275. 

With these conditions, the entropy of the black hole located at r = is finally given by 

S r= o = 7ta/— I r=0 — a 2 , with (6.25) 

c tjk f. (7 — 7) 2 



-I r =o = Q1Q2Q3 [Qe ~ 27m + — Qi \^kj + jlj)(k k + jh) ^ — QjQk 

We finally need to make sure that there are no closed time-like curves in the solution. 
This is done by imposing e~ 4U and e~ 4U r 2 sin 2 6 — uj 1 to be positive everywhere. Given the 
complicated explicit expressions for this quantities, we have not been able to find analytic 
conditions for their positivity. Nonetheless, a careful numerical analysis in different cases 
shows that regularity near each of the centres seems to be sufficient to have the global 
positivity everywhere. 

Mass, Charges, and angular momentum. One can now give the charges associated 
to our solution. The total electromagnetic charges are given by 

g A = -—L- / dv A = -L=(v Asp (0x = 0)-v Av (9 = ir)). (6.26) 

This gives 

P° = --7= ((Qe + Qe) + &%WQ h + Qk)) , go = , (6.27) 

= 77= {h ( d * + dk) - kj(Q k + Q fc )) , Qi = 77= (Qi + Qi) ■ 

In terms of branes, we recall that p° corresponds to the D6 charge, the p*'s to the D4's, 
the (ft's to the D2's and go to the DO. 

It will also be interesting to look at the individual charges of each centre. They 
are easily computed as in (6.26), by picking a cycle surrounding one of the centres. 13 



13 e 



■9- Pa = 572 (K( 6 = 0, S = tt) - w$(0 = tt)) and = ^«(0 S =0)-w$(0 = 0, E = tt)). 
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The electromagnetic charges of the first and the second centre, denoted respectively as 
(Pa>9aa) and (p%,q B \), are given by 







q A i = 


V2 Qi ' 


P°a = 








q B i = 


— Q. 


p° b = 





Pa = -^/f m { k i ~ -rh ' 1 w L ^) Ql 



R 

1 Mk,. 7-7 



R 



QjjQk- (6.28) 



Having computed these charges, one can easily check that, as expected, the entropy of 
each centre (6.25) takes the conventional single-centre expression in terms of local charges 
and the intrinsic angular momenta, i.e. 



S r=0 = ttV - J 4(<?aA,Pa) - « 2 > ( 6 -29) 

where I±{q,p) is the quartic invariant defined in (2.32). The equivalent formula holds for 
the entropy of the second centre. 

The expression of the local charges also allows us to obtain the constraint giving the 
distance R between the two centres as a function of them. Its explicit value is a rather 
complicated function of the charges and the asymptotic moduli. Defining the quantity 

p _ q B iP B pT p { b 2 qAiP A | pT | p a +2 ^ 3Q ^ 

qBi+iq B i+2 q B i+2 q B i+l qAi+iqAi+2 qAi+2 <?Ai+l 

where no sum over % is involved, one can show that 

2(7-7)^ + ^1) =#i, (6.31) 

for all i. This gives that 

D /K qi+l^i+2 — qi+2$i+l (a qo n 

R = — ~n 7 o , (6-32) 

whatever the choice of i is, and with no sum on i. Note that this equation implies 
therefore two extra conditions on the charges of the constituent black holes and the 
asymptotic moduli. Together with the condition that go — 0, the system requires three 
constraints which correspond to the fix value of the phases ctj defining h = 2 Hi in 
this duality frame. It is therefore difficult to interpret the constraint (6.32) as coming 
from the absence of time-like closed curves in this duality frame, and the definition of 
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the correct integrability equation, or bubble equation [5, 31, 37, 38, 39], would require to 
consider the general solution with three arbitrary phases. 

We now want to compute the mass and angular momentum of the solution. In order 
to do this, it is convenient to solve for h in function of the other constants by imposing 
asymptotic flatness, i.e. U — > 

h = 1 , + , T° ~ E kh+ih+2 • (6.33) 

'1*2*3 

i 

For convenience we also define the quantities fa 

ft = kil i+1 l i+2 - m , (6.34) 
such that the asymptotic of the scalar fields reduces to 14 



fi + i 

t % = h+ik+2 j — j2 ■ (6.35) 



In particular, the axions vanish for /, = 0. One computes the ADM mass 

m = ±= (-hi 2 hP° + E m + E 1 + f Z fl+2 qi ) ■ (6 - 36) 

v i i 

This formula reproduce the ADM mass as given in (3.34) for the value of the phase ctj 

e -iai = (i - fi)(i + fi+i)(i + fi+2) ^ 37 ^ 



One important thing to note is that the property that the phases only depend on the 
asymptotic moduli tf Q is due to the simplicity of the specific duality frame we chose; and 
that in general these phases depend non-trivially on the electromagnetic charges, as they 
do for single-centre solutions in particular. Although the mass appears to be linear in 
the electromagnetic charges in this particularly simple duality frame, this is a specific 
property of the latter and the mass formula would be highly non-linear for more general 
charges and asymptotic moduli. 

In fact it is important for the stability that this mass formula does not apply for 
the single-centre solutions with electromagnetic charges q A A,PA an d 1bA,Pb, because one 
would then conclude that the mass of the composite solution is equal to the sum of 



'One has reciprocally f t = e 2<#>i a 4 and l t 2 = ^-SL^ + e t^ - a ~ ) ■ 
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the masses of the constituent black holes, and thus that the solution is only marginally 
stable. Note indeed that one cannot cancel the charges of one centre in (6.28) without 
modifying the electromagnetic charges of the other, such that single-centre solutions with 
the charges (6.28) are not themselves solutions of the system in this particular duality 
frame. One computes instead that the mass formula reproduces the value of the non-BPS 
fake superpotential if and only if 

V l = ~\ O+i + a)fl + (f i+2 + s)f±±), (6.38) 

which is the case when R — > oo for s = m — jhhh- For more general charges, and in 
particular when R is finite and the solution is expected to be stable, the mass formula 
(6.36) corresponds to the 'auxiliary fields' dependent generalised fake superpotential [21] 
for which the auxiliary scalars, say /3„ for a = 1,2, are not extremum of the generalised 
fake superpotential, ^ ^- The auxiliary scalars (3 a parametrize the three 

phases satisfying that (3.34) is real and positive, and (3^ are determined by the con- 
ditions (6.37), and therefore only depend on the asymptotic value of the scalar fields, 
and not the electromagnetic charges. The true 'fake superpotential' (which defines the 
ADM mass of a single-centre non-BPS black hole) is instead defined as W(f3*), for which 
dgff ^ I B-8* = ^" ^ s advocated in [21], W((3*) is bounded from below as a function of 
the moduli, in which case it coincides with the expression determined in [20], whereas it 
develops unstable directions for other values of the auxiliary scalars (3 a . It must therefore 
exist configurations such that W((3 S ) < W(j3*), but this is not true in general. In the 
present situation, the solution will be stable if 

M = W(f3 s ,q A + q B ,p A + p B ) = W(f3 s ,q A ,p A ) + W(f3 s ,q B ,p B ) 

< W((3* A ,q A ,p A ) + W((3* B ,q B ,p B ) = M A + M B . (6.39) 

This is trivially satisfied if both /3 A defines a maximum of W(/3, q A , p A ) and /?* a maximum 
of W(/3, q B ,p B ), but this inequality does not seem to be satisfied for arbitrary charges and 
moduli, as opposed to the BPS case for which the triangle inequality \Z(q A +q B , p A +p B ) \ < 
\Z(q A ,p A ) \ + \Z(q B ,p B )\ holds in general. 

The angular momenta of the solution is given by 



J= -(a + 5) + + ^^)Q 3 Q k 



i-ijikf, , Qi + Qi 

R 

-(a + a)+p^q BA -q AA Pg . (6.40) 



where in the second line we used the definition of the individual charges of each centre. It 
is very interesting to note that, when written in terms of the physical charges, the angular 
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momentum takes the very natural form of the sum of the intrinsic angular momenta plus 
the contribution coming from the electromagnetic interaction between the centres. 
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A Conventions for 50(4, 4) 



In this appendix we give the explicit form of some of the conventions presented in sec- 
tion 3. 

First of all, our representation of 5*0(4, 4) is defined as the one preserving the matrix 
(3.5), it is explicitly given by 
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V o 
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o \ 
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o / 



(A.l) 



The algebra elements appearing in (3.6) 
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V = exp [( A ^E AlB2C:i + aE] exp[UB] 
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The solution of (3.19) in terms of the matrix (3.22) is 



DO : 


(1,1,1,1) 


-N = 


1 — 'j 


= Qi = -P° 








D2 : 


(1,-1,1,1) 


M = 


-Si 


= S2 = S3 = 










(1,1,-1,1) 


M = 


Si = 


— S 2 = S3 = 


1 p2 

2 








(1,1,1,-1) 


M = 


Si = 


S2 = —S3 = 


1 p3 

2 






DA : 


(1,1,-1,-1) 


-N = 


—Si 


= =2 = =3 = 


-Qi = g 2 = 


Q 3 = 


P 




(1,-1,1,-1) 


-N = 


Si = 


— =2 = =3 — 


Qi = -Q 2 = 


Q 3 = 


P 




(1,-1,-1,1) 


-N = 


1=4 = 


=2 = — =3 = 


Qi = Q2 = - 


■Qa = 
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D6 : 


(1,-1,-1,-1) 


M = 


-s, 


— \Qo 








DO : 


(-1,-1,-1,-1) 


-N = 


1 — 'i 


= -Qi = P° 
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— S2 = s 3 = 


2 J 
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— s 2 = s 3 = 


1 p2 
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Si = 


s 2 = —S3 = 


1 p3 
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DA : 


(-1,-1,1,1) 
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Qi = -Q2 = 
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(-1,1,-1,1) 


-N = 




— =2 = =3 = 


-Qi = Q2 = 
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(-1,1,1,-1) 
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=2 = — =3 = 


-Qi = -Q2 


= Q 3 




D6 : 


(-1,1,1,1) 


M = 


-s, 


= -|Qo 









(A.5) 



These definitions deserve some explanations. The four-uplet of numbers give the weight of 
each generator with respect to the four generators if a. In particular, given these numbers, 
one can just read the weight of the generators with respect to h = b A H\. In particular, for 
the BPS orbit, b A = (2, 0, 0, 0) and thus the weight is the first number of the four-uplet. 
In the duality frame that we have chosen in this paper, one can then show that each of 
the generator of the first part of the list, that all have weight one, couple respectively 
to DO, D2, D4, and D6-brane charges. The D-brane generators are then naturally the 
dual one, with opposite weights. From this weights, it is then straightforward to see that 
the almost-BPS orbit (6 A = (4,2,2,2)) corresponds to D6-D4-D2-D0 solutions and the 
composite non-BPS (b A = (0, 2, 2, 2)) to D6-D4-D2-D0 ones. 

Finally, we give here the explicit form of the _K?T-rotation matrix, that takes us from 
one duality frame to another, an element R e KK acts on VV* as 

VV* — >FLVV X R T (A.6) 
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where R is defined as 



R 



R±(a/ , a[) 




R 4 (a' 2 ,a' 3 ) 







) 



(A.7) 



with 



/ cos ccj cos OLj cos ccj sin ctj sin ccj cos atj sin ccj sin atj \ 

— cos ccj sin (X,- cos cos <x,- — sin «j sin atj sin cos ctj 

— sin ccj cos CKj — sin sin ccj cos ai cos ccj cos a« sin «j 

\ sin oij sin «j — sin oij cos «j — cos «j sin «j cos cos a>j / 



(A.8) 



The angles a' A appearing in the R matrix are related to the one of (3.28) through 



References 

[1] A. Strominger, C. Vafa, "Microscopic origin of the Bekenstein-Hawking entropy," 
Phys. Lett. B379 (1996) 99 [hep-th/9601029] . 

[2] R. Emparan and G. T. Horowitz, "Microstates of a neutral black hole in M theory," 
Phys. Rev. Lett. 97 (2006) 141601 [hep-th/0607023]. 

[3] A. Dabholkar, A. Sen and S. P. Trivedi, "Black hole microstates and attractor with- 
out supersymmetry," JHEP 0701 (2007) 096 [hep-th/0611143] . 

[4] D. Astefanesei, K. Goldstein and S. Mahapatra, "Moduli and (un)attractor black 
hole thermodynamics," Gen. Rel. Grav. 40 (2008) 2069 [hep-th/0611140]. 

[5] F. Denef, "Supergravity flows and D-brane stability," JHEP 0008 (2000) 050 
[hep-th/0005049]. 

[6] A. Ceresole, G. Dall'Agata, "Flow equations for non-BPS extremal black holes," 
JHEP 0703 (2007) 110 [hep-th/0702088] . 

[7] G. Lopes Cardoso, A. Ceresole, G. Dall'Agata et al, "First-order flow equations 
for extremal black holes in very special geometry," JHEP 0710 (2007) 063 , 
0706.3373 [hep-th] . 





(A.9) 



46 



[8] A. Ceresole, G. Dall'Agata, S. Ferrara and A. Yeranyan, "First order flows for M = 
2 extremal black holes and duality invariants," Nucl. Phys. B 824 (2010) 239 , 
0908.1110 [hep-th]. 

[9] K. Goldstein and S. Katmadas, "Almost BPS black holes," JHEP 0905 (2009) 058 
, 0812.4183 [hep-th]. 

[10] I. Bena, G. Dall'Agata, S. Giusto, C. Ruef and N. P. Warner, "Non-BPS black rings 
and black holes in Taub-NUT," JHEP 0906 (2009) 015 , 0902.4526 [hep-th] . 

[11] I. Bena, S. Giusto, C. Ruef and N. P. Warner, "Multi-center non-BPS black holes: 
the solution," JHEP 0911 (2009) 032 , 0908.2121 [hep-th]. 

[12] I. Bena, S. Giusto, C. Ruef and N. P. Warner, "Supergravity solutions from floating 
branes," JHEP 1003 (2010) 047 , 0910.1860 [hep-th]. 

[13] S. Giusto, G. Dall'Agata and C. Ruef, "U-duality and non-BPS solutions," 
1012.4803 [hep-th]. 

[14] I. Bena, S. Giusto, C. Ruef, "A black ring with two angular momenta in Taub-NUT," 
1104.0016 [hep-th]. 

[15] K. Goldstein, N. Iizuka, R. P. Jena and S. P. Trivedi, "Non-supersymmetric attrac- 
tors," Phys. Rev. D 72 (2005) 124021 [hep-th/0507096] . 

[16] K. Hotta and T. Kubota, "Exact solutions and the attractor mechanism in non-BPS 
black holes," Prog. Theor. Phys. 118 (2007) 969 , 0707.4554 [hep-th]. 

[17] E. G. Gimon, F. Larsen and J. Simon, "Black holes in supergravity: the non-BPS 
branch," JHEP 0801 (2008) 040 , 0710.4967 [hep-th]. 

[18] E. G. Gimon, F. Larsen and J. Simon, "Constituent Model of Extremal non-BPS 
Black Holes," JHEP 0907 (2009) 052 , 0903.0719 [hep-th]. 

[19] D. Gaiotto, W. W. Li and M. Padi, "Non-supersymmetric attractor flow in symmet- 
ric spaces," JHEP 0712 (2007) 093 , 0710.1638 [hep-th] . 

[20] G. Bossard, Y. Michel and B. Pioline, "Extremal black holes, nilpotent orbits and 
the true fake superpotential," JHEP 1001 (2010) 038, 0908.1742 [hep-th] . 



47 



[21] A. Ceresole, G. Dall'Agata, S. Ferrara and A. Yeranyan, "Universality of the su- 
perpotential for d = 4 extremal black holes," Nucl. Phys. B 832 (2010) 358 , 
0910.2697 [hep-th]. 

[22] J. Perz, P. Smyth, T. Van Riet and B. Vercnocke, "First-order flow equa- 
tions for extremal and non-extremal black holes," JHEP 0903 (2009) 150 , 
0810.1528 [hep-th]. 

[23] S. S. Kim, J. L. Hornlund, J. Palmkvist and A. Virmani, "Extremal solu- 
tions of the S 3 model and nilpotent orbits of G 2{2 )," JHEP 1008 (2010) 072 , 
1004.5242 [hep-th]. 

[24] P. Galli, K. Goldstein, S. Katmadas and J. Perz, "First-order flows and stabil- 
isation equations for non-BPS extremal black holes," JHEP 1106 (2011) 070 , 
1012.4020 [hep-th]. 

[25] P. Galli, T. Ortin, J. Perz and C. S. Shahbazi, "Non-extremal black holes of M = 2, 
d=4 supergravity," 1105.3311 [hep-th]. 

[26] P. Breitenlohner, D. Maison and G. W. Gibbons, "Four-dimensional black holes 
from Kaluza-Klein theories," Commun. Math. Phys. 120 (1988) 295. 

[27] J. L. Hornlund, "On the symmetry orbits of black holes in non-linear sigma models," 
1104.4949 [hep-th]. 

[28] G. Bossard, H. Nicolai and K. S. Stelle, "Universal BPS structure of stationary 
supergravity solutions," JHEP 0907 (2009) 003 0902 . 4438 [hep-th] . 

[29] G. Bossard, H. Nicolai, "Multi-black holes from nilpotent Lie algebra orbits," Gen. 
Rel. Grav. 42 (2010) 509-537. , 0906 . 1987 [hep-th]. 

[30] G. Bossard, "1/8 BPS black hole composites," 1001.3157 [hep-th] . 

[31] B. Bates and F. Denef, "Exact solutions for supersymmetric stationary black hole 
composites," [hep-th/0304094]. 

[32] D. Collingwood and W. McGovern, "Nilpotent orbits in semisimple Lie algebras" 
Van Nostrand Reinhold Mathematics Series, New York, 1993. 

[33] J. Sekiguchi, "Remarks on real nilpotent orbits of a symmetric pair," 
J. Math. Soc. Japan 39 (1987) 127138. 



18 



[34] D. Z. Bokovic, N. Lemire, J. Sekiguchi, "The closure ordering of adjoint nilpotent 
orbits in So(p,q), Tohoku Mat. J. 53 (2001) 395. 

[35] F. Denef and G. W. Moore, "Split states, entropy enigmas, holes and halos," 
[hep-th/0702146]. 

[36] K. Behrndt, R. Kallosh, J. Rahmfeld et al, "STU black holes and string triality," 
Phys. Rev. D54 (1996) 6293. [hep-th/9608059] . 

[37] I. Bena, N. P. Warner, "Bubbling supertubes and foaming black holes," Phys. Rev. 
D74 (2006) 066001. [hep-th/0505166] . 

[38] P. Berglund, E. G. Gimon, T. S. Levi, "Supergravity microstates for BPS black holes 
and black rings," JHEP 0606 (2006) 007. [hep-th/0505167] . 

[39] A. Saxena, G. Potvin, S. Giusto, A. W. Peet, "Smooth geometries with four charges 
in four dimensions," JHEP 0604 (2006) 010. [hep-th/0509214]. 



19 



